Criteria for quantum coherent transfer of excitons between chromophores in a polar 

solvent 
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We show that the quantum decoherence of Forster resonant energy transfer between two opticaUy 
active molecules can be described by a spin-boson model. This allows us to give quantitative 
criteria, in terms of experimentally measurable system parameters, that are necessary for coherent 
Bloch oscillations of excitons between the chromophores. Experimental tests of our results should 
be possible with Flourescent Resonant Energy Transfer (FRET) spectroscopy. Although we focus 
on the case of protein-pigment complexes our results are also relevant to quantum dots and organic 
molecules in a dielectric medium. 



Decoherence is the process whereby quantum inter- 
ference effects are "washed out" by the interaction of a 
quantum system with its environment. It has been sug- 
gested that decoherence is responsible for the crossover 
from quantum to classical behavior Decoherence 
places limits on the possibility of quantum computation 

• The challenge of building a quantum computer and 
the potential of biomimetics raises the possibility of 
exploiting the self assembly of complex biomolecular sys- 
tems, such as light harvesting photosynthetic protein- 
pigment complexes However, this also raises pro- 
found questions, which are of interest in their own right, 
about what role quantum effects play in biomolecular 
functionality. One model for describing quantum deco- 
herence is the spin-boson model 0, |^ . It describes quan- 
tum tunneling between two quantum states that are cou- 
pled to a dissipative bath which is modelled by a set of 
harmonic oscillators (see the Hamiltonian in eqn. H1U|) 
below). This model has been used to describe systems 
ranging from Josephson junction qubits to electron 
transfer in biomolecules In this Letter, we show 

how the spin-boson model can also be used to describe 
the transfer of excitons between two chromophores by 
the mechanism of Forster resonant energy transfer in a 
dielectric medium. Established results for the spin boson 
model are then used to give stringent criteria for quantum 
coherent exciton transfer. Our results are also relevant to 
quantum dots (see for example ^^'^ small molecules 
in a dielectric medium fllj . 

Model for interaction of the individual chromophores 
with the solvent. In can be shown [l^ that the coupling 
of the electronic excitations in a chromophore to its en- 
vironment may be modelled by an independent boson 
model lial of the form 
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FIG. 1; Two chromophores with energy gaps ei and £2 are 
coupled with an interaction energy A due to the Forster 
dipole-dipole interaction. In the simplest model, the chro- 
mophores are centred inside spherical proteins with static di- 
electric constant ep, and surrounded by a polar solvent with 
frequency dependent dielectric constant e{u)). 



experienced by the chromophore dipole due to the "cage" 
of polarised solvent and protein molecules around it. The 
Ca are the couplings of the excitation to each mode. The 
coupling to the environment and quantum dynamics is 
completely specified by the spectral density. 
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In the simplest picture of protein-pigment complexes, 
the chromophore can be treated as a point dipole in- 
side a uniform, spherical protein 0,0] surrounded by a 
uniform polar solvent |l4| . Typical dielectric relaxation 
times of proteins ^3] are significantly longer than the 
other time relevant scales (except fluorescence lifetimes, 
see Table nj and so we consider only the static dielectric 
constant for the protein. The spectral density ^3 is then 
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Here the chromophore is treated as a two level system 
with energy gap e between the ground and excited state, 
and A^ is the difference between the dipole moment of 
the chromophore in the ground and excited states. R = 
^^Ca{aa + ajj) is the quantised reaction field 0,0] 



where h is the radius of the protein containing the chro- 
mophore, e(w) is the complex dielectric function of the 
solvent and is the static dielectric constant of the pro- 
tein. For the case of a Debye solvent [l^ . 
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where eg and are the static and high frequency fre- 
quency dielectric constants of the solvent respectively, 
and te — ^^^t^td where tjj is the Debye relaxation 
time of the solvent. For water at room temperature, these 
parameters are Es = 78.3, eoo = 4.21 and Tjy = 8.2ps 
[l9j while for THF (tetrahydrofuran) they are Cs = 8.08, 
e.i = 2.18 and Tij = 3ps Typical protein static dielec- 
tric constants Cp are between 4-40 dep ending on which 
part of the protein is of interest [l^ [I^, ■ therefore 
takes values between 0.5 — 2.5ps. 

Model for FRET in the presence of a solvent. We 
now consider the case of two biomolecules coupled by 
the Forster interaction . This is a dipole-dipole inter- 
action which produces a non-radiative transfer of an ex- 
citation between two chromophores (see Figure^. This 
interaction is the basis for energy transportin photosyn- 
thetic light harvesting complexes and fluorescent reso- 
nant energy transfer (FRET) spectroscopy. We may 
write the total Hamiltonian as the sum of two spin-boson 
Hamiltonians for each chromophore |l2|: 
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where R, = J2a^'haiai,a 



{i = 1,2) is the quan- 



tised reaction field operator for molecule i, and Bi — 
J2a ^i,aol a^i,a IS the energy stored in the solvent cage of 
molecule i. For molecules sufficiently far apart (> 20 A), 
the interaction is given by 
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where = {e\fl\g)^, i = 1,2, is the transition dipole 
moments of the chromophores (distinct from the change 
in dipole moment of the molecule during the transition, 
A/x;), n is the refractive index of the solvent, R is the 
separation of the molecules and k is related to the relative 
orientation of the two dipoles . 

For future convenience, in matrix notation the Hamil- 
tonian (0) 



H = 
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(7) 



/e+ + V+ 



V 





e-+V- 
A 





A 





where 
V± = AfiiRi ± Afi2R2 



Ami C„(a„ + al) ± A^a + ^p) 

a /3 
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TABLE I: Typical energy scales E and the corresponding time 
scales T — h/E relevant to Forster Resonant Energy Trans- 
fer in several systems. A is the Forster coupling between the 
two chromophores, e is the difference in energy of their first 
excited states. lJc is the high frequency cut-off of the spec- 
tral density J{u!) in the relevant spin-boson model. Typically, 
I/lOc = -^^^^'''D, where td is the Debye relaxation time of 
the solvent and es and ei are the solvent's static and high fre- 
quency dielectric constants respectively. LH-I and LH-II are 
the photosynthetic light harvesting complexes, and "BChl" is 
the bacteria chlorophyll molecule. We observe that the flu- 
orescence lifetime Tj-^^j is much longer than the other time 
scales, suggesting that all other processes of interest occur 
before radiative decay. We note also that both A and e span 
two orders of magnitude, so we might expect very different 
behaviour for BChl's in LH-II than for typical FRET spec- 
troscopy, such as between green and red fluorescent proteins. 



System 


E (meV) 


r (ps) 


Ref 


300K 




25 


0.16 




H2O J{uj) cut-off 




2-8 


0.5-2 


fl2] 


THF J(w) cut-off 




2-4 


1-2.5 


[20] 


Radiative lifetime 




4 X 10"* 


10* 


[24] 


Protein relaxation time 


h/Tp 


0.025 


162 


fl7] 


Typical FRET 


A 


0.2-2 


2-20 


[25] 


(Green — > red) 


e 


500 


8 X 10"=* 




BChl in LH-II 


A 


46-100 


0.04 - 0.08 


[26] 




e 







f26j 


LHII -> LHII 


A 


0.3 


13 






e 







m 


LHII LHI 


A 


0.6 


7 


[26] 




e 


6.5 


0.6 


[26] 



and e± = ei ± £2- 

Mapping to the spin boson model. The number of ex- 
citations in the system is related to iV = 5(0-^ + a"^ + 2) 
which we note commutes with the Hamiltonian: [H, N] — 
0. Hence the total number of excitations is a constant 
of the motion. Note we are assuming that the fluo- 
rescence lifetime of the chromophores is much longer 
than the other time scales of the system described by 
H and so do not need to include radiative decay in 
H (typically 1 — 10ns, Table QJ. If we consider only 
singly excited systems then TVlvP) = |^), and we can 
project onto the corresponding two dimensional subspace 
15)2 , 15)1 |e)2}, where \g)^ , |e)^ represents the 
ground and excited states respectively of chromophore 
i = 1,2. We note that this is a decoherence-free subspace 
with respect to decoherence due to the environment. 
We can thus restrict our Hamiltonian to the central 2x2 
submatrix of Q), which can be written in terms of new 
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Pauli sigma matrices as the spin-boson model 



(9) 



where V = represents the interaction with the en- 
vironment, and e = e_ ~ ^ ^2 is the difference in 
energy for the excitation on the different chromophores. 
We assume the bath modes coupled to each chroniophore 
are independent, i.e., [ai,Q,,a2^] = 0. This can be justi- 
fied if the molecules are sufficiently far apart and their 
cavities can be treated independently |27l |. The environ- 
ment can then again be modelled as a set of independent 
harmonic oscillators in the standard form of the spin- 
boson model: 



H 



1 



(10) 



where the Oq now include both sets of independent har- 
monic oscillators. 

To complete the description, we must specify the new 
spectral density J(w) which now describe the environ- 
ments around both molecules jointly. This may be ob- 
tained as for the single molecule case J^j by following the 
ansatz of Caldeira and Leggett by treating the fluc- 
tuations in the environment via the correlation function 

{vmrn-. 

{vmm = 



A^iA^ia (i?iWi?2(0)) + {R2{t)Ri{0) 



Provided that the chromophores are sufficiently far 
apart that their cages are uncorrelated, (i?i(t)i?2(0)} = 
(i?2(i)^i(0)) = 0. J{uj) is then given by 



J(cj) = Im J dte^' 



[i{v{t)vme{t)] (11) 



Jl{uj)+ J2{lo), 



(12) 



i.e., the new spectral density is simply the sum of the 
appropriate spectral densities for the individual chro- 
mophores. 

Modelling of FRET systems by spin-boson like mod- 
els has been considered previously [23, ll^, |33 , but only 
with perturbation theory and Fermi's golden rule. Here 
we have presented a specific microscopic derivation of 
the effect of the environment on transfer, and given an 
exphcit form for the spectral density J[ijj) that can be 
determined from experiments on single chromophores. 

We note that we have Ohmic dissipation, i.e., J(w) = 
aoj for Lo < LOc — H/te with the dimensionless coupling 
constant 
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FIG. 2: Dependence on the environmental coupling a of the 
critical temperature T* for the cross-over from coherent to 
incoherent behaviour. Above the sharp cut-off at a = 1/2, 
coherent oscillations do no occur for any temperature. Re- 
produced from ref. 's]. 



If e, A <^ LUc, for Ohmic dissipation a is a critical pa- 
rameter for determining the quantum dynamics |5|, |a, l3lj 
(see Table . For typical free chromophores in water 
at room temperature, a of the order 0.1 — 1 [l^ . which 
represents strong coupling to the environment (in com- 
parison, a is orders of magnitude smaller for Josephson 
Junction qubits 7].) However, the protein environment 
pushes the solvent away from the chromophore, and can 
make a much less than one. 

Criteria for quantum coherent RET. The location of 
the excitation is given by P{t) — {az{t)). Suppose that 
the excitation is initially localised on one chromophore 
(the "donor"), corresponding to P(0) ~ +1, with the 
other chromophore (the "acceptor" ) in the ground state. 
The Forster coupling will cause transfer of the excita- 
tion between the chromophores. We now establish the 
conditions of validity of Forster's equation for FRET ef- 
ficiency in terms of the convolution of the ab sorp tion and 
emission spectras of the two chromophores [j^^Evj such 
as is widely used in "spectroscopic ruler" applications 
[3^ in molecular biophysics. These results are based 
on two assumptions (i) second order perturbation the- 
ory in A (the Fermi Golden Rule), which assumes that 
A is small compared to the other energy scales of the 
system, and (ii) that there is no back transfer, i.e., inco- 
herent transfer of the excitation. In typical FRET sys- 
tems, such as between red and green fiuorescent proteins, 
A is typically less than ImeV, while e ~ 500meV and 
hwc ^ 8meV (see Table IJ) which justifies Forster's as- 
sumption (i). However, assumption (ii) is only justified 
by Leggett et al.'s nontrivial derivation: provided that 
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A < Hu>c, then for a sufRciently large difference in energy 
levels (~ e > = A(A/a;c)"/^^""\ where A^ is the os- 
cillation frequency renormalised by interaction with the 
bath environment) or for sufficiently high temperature 
(T > T* - Ar/a, see Table U and Figure 13 the transfer 
is always incoherent |^. We have therefore justified the 
use of Forster's equation for typical FRET systems. 

However, when e is small and < a < 1/2 and 
T < T*{a) (Table ini) coherent oscillations may occur 
d IHJ and the Fermi Golden Rule derivation for FRET 
no longer applies. Further, from Tabled we see that A 
may be comparable to or greater than the reorganisation 
energy hojcCe [^^^. Therefore, second-order perturba- 
tion theory in A may be hard to justify. From Table 
m we see that A takes on a wide range of values, and 
it is conceivable that at low temperatures and between 
chromophores with very close (or identical) energy lev- 
els that the assumption of incoherent transfer, and hence 
the Fermi Golden Rule result, may break down. 

Experimental tests. A possible way to observe the co- 
herent oscillations of excitons between chromophores is 
to use identical chromophores which are at an angle to 
each other so that their dipole moments are not paral- 
lel. Then from a bulk sample the oscillations should be 
present in the time dependence of both the fluorescent 
anisotropy and the flourescent noise js^ . This effect has 
recently been reported for pairs of chromophores (e.g., 
anthracene dimers) that are covalently bonded and in a 
solvent at room temperature [s^lls^ l. However, although 
in all cases the conditions of the experiment are in the 
regime A < Wc they do not satisfy two of the necessary 
conditions for coherent oscillations, T < T* ~ A/a and 
a < 1/2. For example, reference [221 concerns a molecule 
DTA in THF (tetrahydrofuran) solvent at room tempera- 
ture. The measured oscillation period is 1 picosecond and 
the damping time is also 1 picosecond. For this solvent 
the cut off frequency is — 4 meV (see Tabled), and we 
expect A < hiUc and so the results for Ohmic dissipation 
should be relevant at the qualitative level. To determine 
a, we can estimate the reorganisation energy, Er ~ ahiUc, 
from the width of the absorption and fluorescent spectra 
or from the Stokes shift Roughly, - 100 meV, 

and so a ~ 25 and no coherent oscillations should occur. 
Furthermore, Monte Carlo simulations that take into ac- 
count a broad range of A /ojc values do not give coherent 
oscillations for these kind of parameter values (compare 
Fig. 13 in ref [Sfl and Fig. 7 in ref jH). 

Application to photosynthesis. One system of particu- 
lar interest is the transfer of excitations between bacterial 
light harvesting (LH) complexes (I and II) in photosyn- 
thetic units. In a typical process, an LH-II ring of chloro- 
phyll chromophores absorbs a photon. The excitation 
may then transferred to other LH-II rings before reach- 
ing the LH-I ring where it is sent to the reaction centre to 
be converted to chemical energy. As the chromophores 
within the ring are identical, e « while A is between 



TABLE II: Beliaviom- of P{t) = {cr^{t)), which gives the lo- 
cation of the excitation as a function of time t, for e = and 
A <^ hiUc, where A is the Forster coupling strength and lJc is 
the high frequency cut-off of the spectral density J{u>). "loc" 
refers to localisation, "coh" to damped coherent oscillations 
and "inc" to incoherent behaviour i.e., exponential decay. T 
is the temperature of the system and a is the dimensionless 
coupling constant of the chromophores to the environment as 
defined in Equation 1131 . r refers to the relaxation rate in the 
expression P{t) — exp(— t/r). The analytic form of P{t) is 
given where known, and is generally valid only for timescales 
longer than I/lOc. Ar = A{A/lj^)°/'-^-"'> 



a> 1 
< 1/2 

< 1/2 
> 1/2 



T 



T = 
T=0 

T >T* 
T > 



Key 



loc 
coh 

inc 
inc 



P{t) 



P{t) = 1, all t 
P{t) ^ exp [-2t 

X cos sm 

-1 _ A'^ r(a) r -n-fcT l 

2 r(Q+l/2) [ hLJa \ 



— sin^ , 



1-a 



Ref 



\5] 
[31] 



46 — lOOmeV. Here A ^ e, but also A > uic, and so 
the results of are not applicable. We expect coherent 
transfer of the excitation around the ring [s^ . Experi- 
mental studies suggest that the excitation is indeed de- 
localised ll^. For transfer between identical LH-II rings, 
again e = but A <^ hujc (see Table IJ). In this case, 
T > T*{a) (Figure and the transfer will be incoher- 
ent. Finally, for transfer on an excitation from an LH-II 
to an LH-I ring, e 3> A and we expect incoherent trans- 
fer. These behaviours play an important functional role 
- the delocalisation of the excitation prevents radiative 
loss of the excitation [2^ , while the incoherent inter-ring 
transport establishes a one-way flow of energy towards 
the reaction centre. 

In conclusion, we have shown that the decoherence of 
two chromophores coupled by the Forster interaction in 
the presence of a solvent can be described by the spin- 
boson model. We find that while the standard FRET ef- 
ficiency formulas do hold for most typical systems, they 
are based on nontrivial assumptions which may not be 
true in general. We give the quantitative conditions nec- 
essary for coherent Bloch oscillations of excitons between 
the chromophores and suggest how this could be tested 
experimentally. Finally, we have used these models to de- 
scribe the transfer of excitations between light harvesting 
complexes in bacteria and their relevance to the system's 
biological functionality. 
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